Towards a relevant set of state variables to describe static granular packings 
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We analyze, experimentally and numerically, the steady states, obtained by tapping, of a 2D 
granular layer. Contrary to the usual assumption, we show that the reversible (steady state branch) 
of the density-acceleration curve is nonmonotonous. Accordingly, steady states with the same 
mean volume can be reached by tapping the system with very different intensities. Simulations 
of dissipative frictional disks show that equal volume steady states have different values of the 
force moment tensor. Additionally, we find that steady states of equal stress can be obtained by 
changing the duration of the taps; however, these states present distinct mean volumes. These 
results confirm previous speculations that the volume and the force moment tensor are both needed 
to describe univocally equilibrium states in static granular assemblies. 
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PACS numbers: 

Introduction: Finding out the appropriate set of 
macroscopic variables that characterizes the equilibrium 
state of a given system is the first step in any thermody- 
namic study [l[ . What these variables are for a sample of 
sand in a box is still under discussion Q. Twenty years 
ago, Edwards and Oakshott[l] put forward the idea that 
the number of grains N and the volume V are the basic 
state variables that suffice to characterize a static gran- 
ular sample in equilibrium. The NV granular ensemble 
was introduced as a collection of microstates, where the 
sample is in mechanical equilibrium, compatible with N 
and V. 

Experimentally, such equilibrium states are commonly 
obtained by tapping the system at different dimension- 
less accelerations, T, which is the standard parameter 
used to quantify the intensity of the external excitation 
Hence, under the action of repeated excitations at 
constant T, the system properties reach a steady state 
where all observables fluctuate around well defined mean 
values. Therefore, tapping provides the energy necessary 
to allow the system to explore different configurations 
over which macroscopic observables can be averaged. In 
the steady state, the mean volume V, or, more commonly, 
the mean packing fraction, <p (defined as the percentage of 
the space occupied by the particles), can be measured. It 
is now well established 4-8] that different steady states 
can be reached reversibly by varying V. These experi- 
ments also suggest a monotonic relation between T and 
<f). Nevertheless, (j> does not only depend on T. In har- 
monic pulses, if r is kept constant while the duration 
and the amplitude of the pulse are simultaneously mod- 
ified, the system evolves to different steady states. This 
fact was used to show that equivalent equilibrium states 
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[8J, |9( can be generated by using different combinations of 
pulse amplitude and duration. On the other hand, new 
theoretical works [Tol - fTB ] suggest that the force moment 
tensor, E, (E = Va, where a is the stress tensor) must be 
added to the set of extensive macroscopic variables (i.e., 
an NVE ensemble) to describe adequately a packing of 
grains. 

Recent numerical simulations fltjl \v\ have predicted 
also the existence of a nonmonotonic relation between T 
and <p, if sufficiently large values of T are considered. A 
hint as to this behavior can aheadybe seen in some pack- 
ings obtained in the laboratory [18|; although this has 
passed unnoticed by those studying equilibrium proper- 
ties. These results imply that equilibrium states of equal 

V (or <fi) could be obtained not only by changing the pulse 
duration, as in Ref. Q , but also by simply using different 
values of the excitation parameter T. In what follows, we 
will show that these equal volume equilibrium states are 
experimentally accessible. We will then use numerical 
simulations to show that these states are, nevertheless, 
different. We have found that two equilibrium states ob- 
tained with different excitations, even if they have same 

V and N, may present distinct values of E. This result 
suggests the incomplete picture upon which research has 
been based for over two decades when the properties of 
static granular matter have been investigated; and takes 
us one step closer to completing the definition of the state 
variables that are necessary to describe a static packing 
of grains. 

Experimental results: A quasi 2D Plexiglass cell 
(width: 28 mm, height: 150 mm) was used to study the 
packing dynamics. The cell was filled with 1000 alu- 
mina oxide (AI2O3) spheres of diameter 1.000 ± 0.005 
mm. The separation between the Plexiglass sheets was 
10% larger than the bead diameter in order to mini- 
mize the particle- wall friction and prevent arching in the 
transversal direction. The system was tapped with an 
electromagnetic shaker (Tiravib 52100) with a series of 
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FIG. 1: (Color online) (a) Evolution of <f> towards the steady 
state at r = 7.0 starting from the steady state correspond- 
ing to r = 1.5. (b) The evolution when switching back to 
r = 1.5 (Arrows indicate the tap where F is switched), 
(c) Steady state packing fraction, <j>, as a function of F (green 
down triangles: steady state obtained in Exp A through tap- 
ping at constant F from an initial ordered state; orange up 
triangles: steady state in Exp B obtained by an annealing 
protocol) 0, [2lJ . (d) Steady state relative packing fraction 
4>/4>ma,x as a function of tapping intensity F for the 3D cell 
(0max corresponds to the value of <j> obtained at the lowest F 
studied). Error bars represent the standard error of the mean 
estimated from the standard deviation. 



harmonic pulses of variable amplitude and constant fre- 
quency v = 30 Hz every three seconds. The tapping in- 
tensity was measured with a piezoelectric accelerometer 
attached to the base of the cell. Although in recent years 
alternative parameters have been proposed to quantify 
the intensity of the taps [l9|, US], we employ the usual 
nondimensional peak acceleration T = a pea k/g (where 
g is the acceleration of gravity) . High resolution digital 
images of the packing were taken after each tap. The cen- 
ter of each sphere was detected with an error of less than 
2%. The packing fraction of each image was obtained 
by considering each grain as a disk of the corresponding 
effective diameter and then calculating the percentage of 
the area covered by the disks in a rectangular area 10% 
smaller than the size of the packing. Due to the space 
left between the Plexiglass sheets, small overlaps can be 
observed in the 2D projection taken by the photographs. 
This effect leads to estimations of 4> above the 2D hexag- 
onal close packing value for the densest configurations 
studied. After a relative short number of taps, the sys- 
tem reaches a steady state characterized by a plateau in 
curves of cj) vs. the number of taps. In contrast with 
the three dimensional case, any stationary steady state 



is obtained after a few taps. We checked that the mean 
value and standard deviation of the packing fraction dif- 
fer less than 0.1% if 200 or 2000 taps of equilibration 
are applied. As an example, in Fig. 1(a), we show that 
the steady state corresponding to T = 7.0 can be reached 
in a couple taps, even when the initial packing is in an 
ordered configuration. Equilibrating the system back at 
r = 1.5 takes somewhat longer (about 100 taps) [see Fig. 
1(b)] [2l|. Therefore, the sample initially prepared in a 
highly ordered configuration was tapped 500 times for 
any given T before taking averages over 100 taps, and 10 
independent runs were averaged. Alternatively, an an- 
nealing protocol was used in which T was decreased in 
discrete steps from the highest values and tapped 500 
times at each T value without emptying the cell. 

In Fig. 1(c), we plot <f> in the steady state as a function 
of r for different independent repetitions of the experi- 
ment. It can be seen that the same results are obtained 
by equilibrating from an initial ordered structure (Exp. 
A) or by following an annealing path (Exp. B). Although 
the fluctuations of cf> are large [see Fig, 1(a)], its mean 
value is well defined with a small confidence interval [see 
error bars in Fig. 1(c)]. For low excitations, <fi decreases 
as r is increased, in agreement with previous results re- 
ported by several groups. However, beyond a certain 
value r m j n , the packing fraction grows [21| . The same 
trend is observed if the tap frequency v is changed (Exp. 
C) [see Fig 2(a)]. An explanation for this behavior based 
on the formation of arches has been given in jl6| . 

In order to assess whether the presence of a mini- 
mum in the T-cj) curves is induced by the highly ordered 
crystal-like structures present in the quasi 2D cell, we 
repeated the experiment with a 3D cell. This 3D setup 
consists of a cell of the same height and width (as the 
quasi 2D cell) but 6 mm thick. In this case, the granular 
sample was made of polydisperse glass beads 1.0 ± 0.2 
mm in diameter. For the sake of simplicity, in the 3D 
cell, the relative packing fraction was estimated from the 
height of the granular layer. The stationary regime in the 
3D setup was obtained after 2 x 10 4 taps. In Fig. 1(d) 
we show the steady state packing fraction as a function 
of r for this three-dimensional setup. Again, the same 
trend as in the quasi 2D experiment is observed. This 
nonmonotonic dependence challenges the idea that the 
volume defines the equilibrium state, unless these equal 
volume states are proved to be equivalent [24j |. 

Numerical evidences: In 2005, Edwards suggested that 
the stress tensor should be included in the description of 
the equilibrium state of a static granular sample along 
with N and V in the so-called full canonical ensemble 
[Til ] . There is some consensus now that the 

force moment tensor, E = V a , allows for a microcanon- 
ical description [141 ] . However, no experiments or simu- 
lations have shown that this variable is necessary to dis- 
tinguish between different states of thermodynamic equi- 
librium. 

Since it is difficult to measure the force moment tensor 
in our experimental setup, we use a realistic computer 
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FIG. 2: (Color online). Results from the MD simulations of the soft disk model for two different frequencies of the tapping 
pulse [blue squares: v\ = Q.5(g/d) 1 ^ 2 , and red circles: U2 = 0.25((?/e() 1/ ' 2 ]. (a) Steady state packing fraction, <f>, as a function of 
the drive F. In order to compare with the quasi 2D experiments (green down triangles, orange up triangles and black stars) 
the vertical and horizontal axes have been scaled with the characteristic values </> m i n and F m i n , respectively, (b) Fluctuations, 
o qy , of the packing fraction in the steady state (measured by the standard deviation) as a function of </)■ The arrows indicate 
the direction of increasing tapping intensity. Error bars correspond to the standard error, (c) The trace, Tr(E), of the force 
moment tensor in the steady state as a function of r/F m i n . 



model of the quasi 2D setup [H, [22[ . We used a veloc- 
ity Verlet algorithm to integrate the Newton equations 
for 512 monosized disks in a rectangular box. The disk- 
disk and disk-wall contact interaction comprises a linear 
spring-dashpot in the normal direction and a tangential 
friction force that implements the Coulomb criterion to 
switch between dynamic and static friction. Units are 
reduced with the diameter of the disks, d, the disk mass, 
m, and the acceleration of gravity, g. Details on the 
force equations and the interaction parameters can be 
found elsewhere [lH HH . Tapping is simulated by mov- 
ing the confining box in the vertical direction following 
a half sine wave trajectory. The intensity of the exci- 
tation is controlled either through the amplitude, A, or 
the frequency, is, of the sinusoidal trajectory; and it is 
characterized by the parameter T — A(2nv) 2 / g. We con- 
sider the system has reached the steady state whenever 
a plateau in <fi is observed (with no visible trend in a plot 
of <fi vs \og(taps)) . Averages were taken over 400 taps in 
the steady state and over 20 independent simulations for 
each value of A and v. 

In Fig. 2(a), we show </> as a function of T for this 
model. Different values of T were obtained by varying 
both A and v. As in the experiment, states of equal <f> 
are generated at both sides of a minimum packing frac- 
tion, (f> m - m . States of equal volume at each side of (j>min 
display slightly different volume fluctuations [Fig. 2(b)] 
suggesting that these states are not the same. In con- 
trast to the packing fraction, the trace of E [23[ gently 
grows as T is increased [see Fig. 2(c)]. There is neither a 
minimum nor a maximum in the T-T, curves. Therefore, 
states of equal </> at each side of (/> m in present distinct 
£. Hence, states of equal volume at each side of r m i n 
are, in fact, different. This finding proves that the vol- 
ume alone cannot characterize the equilibrium states of 
a static granular sample and suggests that the force mo- 
ment tensor is a convenient extra state variable. 



In Fig. 3(a), the loci defined by the equilibrium 
states visited during the numerical tapping experiments 
are plotted in a hypothetical </>-E thermodynamic phase 
space. It is clear that states of equal E can present dif- 
ferent </> if prepared at different v. Thus, E cannot define 
the equilibrium state by itself. Both V and E have to 
be specified to fully identify a given state. This result is 
confirmed by the fluctuations of the force moment ten- 
sor, <7£, [see Fig. 3(b)]. The fluctuations as for a given 
E are not unique: they depend on the frequency v of the 
tapping. However, for the states obtained with different 
excitations that display equal </> and E (note that such 
states would actually be a unique state), we find that 
ct^ and <T£ display also very similar values. Hence, the 
coincidence of the mean and fluctuations of the state vari- 
ables suggests that no other extensive parameter would 
be necessary to describe such equilibrium states. 

It is interesting to mention that the definition of an 
entropy S(N, V, E) is particularly convenient since E is 
not bounded as V is. This implies that S can be a mono- 
tonic increasing function of E, which permits a well be- 
haved Legendre transformation to a canonical ensemble 
Hence, one can define a non- negative temperature- 
like quantity: the angoricity [Hldil — the proposed name 
for the inverse of the conjugate variable to the force mo- 
ment tensor. Let us finally recall that the temperature- 
like quantity associated with the volume (so-called com- 
pactivity) can present negative values due to population 
inversion [25j ]. 

Conclusions: We have shown that steady states of 
static granular packings obtained by tapping the system 
with different pulse strength and duration can present 
different mean force moment tensor even if they cor- 
respond to the same mean volume. Moreover, steady 
states that present the same mean force moment ten- 
sor are distinguishable by their mean volume. To our 
knowledge, this is the first experimental/numerical evi- 
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dence that both extensive variables should be included in 
an entropic formulation of the thermodynamics of such 
steady states. 
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FIG. 3: (Color online). The loci of the generated equilib- 
rium states in the </>-E phase diagram, (a) The stationary 
packing fraction as a function of the trace Tr(E) for v\ (blue 
squares) and (red circles), (b) The fluctuations, or,.^), 
versus Tr(E). Error bars as in Fig. 2. The shaded areas 
indicate the position of <^ m in- 
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